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Abstract 

We compute the spectrum and the eigenstates of the open XXX model with non- 
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nal, the usual coordinate Bethe ansatz does not work anymore, and we use a general- 
ization of it to solve the problem. 
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1 Introduction 



The XXX model with periodic boundary conditions [T] is one of the most studied model 
in the realm of integrable systems. It was solved by Bethe [2], using what is now called 
the coordinate Bethe ansatz, and since then, many papers appeared on the subject. As for 
most of integrable models, the situation drastically changes when one considers the open 
case, i.e. when the boundary conditions are not periodic anymore. From the pioneer works 
of Cherednik [3] and Sklyanin it is known that integrability of the model is preserved 
when the boundary conditions are coded by two independent matrices obeying the so-called 
reflection equation. Classification of such matrices amounts to classify integrable boundary 
conditions. However, although the models are known to be integrable, there is no general 
answer to get the Hamiltonian eigenstates for generic integrable boundaries (even when the 
periodic ones are known). 

For the XXX model and its su{N) generalization, such classification of boundaries has 
been done in [5H7]. However, as mentioned, the models are still not solved for generic 
boundary matrices. The situation where the two matrices are diagonal is well-understood: 
the spectrum can be computed using analytical Bethe ansatz [8]-[T0]. the eigenfunctions can 
be obtained from coordinate [TTHT3] or algebraic |1| [H] - [T6] Bethe ansat^, and correlation 
functions are well studied, see e.g. [T7HTU] . The case of simultaneously diagonalizable ma- 
trices is done in [2U] (see also [T5]). 

The spectrum for some non diagonal cases, where the boundary parameters (entering the 
boundary matrices) obey relations, could be deduced from studies on XXZ models, as dealt in 
e.g. [2T1[22]. Let us also mention [23] where numerical methods are used to get the spectrum 
and [2^ where a deformation of Onsager algebra is studied to compute the spectrum for 
general boundary matrices. The first attempt to compute the spectrum and eigenfunctions 
for XXX model with boundary matrices that are not simultaneously diagonalizable can be 
found in [25]. There, the algebraic Bethe ansatz is used but one must restrict itself to one 
triangular boundary and one diagonal. Moreover, the su{2) invariance of the i?-matrix is 
needed, so that the treatment is specific to the XXX model. The aim of this article is to 
tackle the same question but using a generalized coordinate Bethe ansatz (gCBA). Indeed, 
recently, in the case of XXZ model such problem has been successfully solved using such 
a method [26l[27]. It leads to a simpler presentation of the gCBA, that we present from 
the integrable models view point (while the presentation in [261127] was more ASEP like). 
We hope that this presentation permits to emphasize the novelties of the gCBA without 
the technical difficulties encountered in the ASEP model. The gCBA allows us to construct 
the eigenstates and compute the spectrum for one diagonal and one non-diagonal boundary 
matrix. The solution is equivalent to the g — 1 limit of the XXZ model [211 [221 [27] 
the one found by algebraic Bethe ansatz [25]. The constraints on the boundary matrices for 
XXZ model found previously [2Tl[22l[27] being now replaced by the triangular form of the 
left boundary matrix and the diagonal form of the right one. 

To be more specific, we will construct the eigenfunctions and compute the spectrum and 

^Strictly speaking, references [9l[T4] deal with the XXZ model. One has to perform a q ^ 1 limit to get 
the XXX model. 
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Bethe equations for the following Hamiltonian: 

L-l L-1 

H = Bt + Hbulk + BJ^ where H^uik = h^j+i = (Pi,£+i - l) (1.1) 

e=i e=i 

with the following boundary matrices 

Let us stress that, generically, both boundary matrices are not simultaneously diagonalizable 
(or even not diagonalizable at all), since they do not commute and that the total spin is not 
a good quantum number, since it does not commute with the Hamiltonian (because of the 
triangular boundary matrix). Of course, the method also applies when <H- B~ or when 
one considers a lower triangular matrix instead of a upper triangular one. Moreover, one 
can also conjugate both B~^ and B~ by the same matrix: the spectrum is the same, and the 
eigenstates are constructed in an obvious way. 

The plan of the paper is as follows. Notations are detailed in section O The construction 
of the eigenstates is done in section [3], and section H] is devoted to the proof of the main 
property of the article. We conclude in section [51 



2 XXX model with boundaries 

The model we consider has Hamiltonian (11.11) . It acts on a spin chain with L sites, with a 
spin I (a space) on each site of the chain. Hence, the Hilbert space is n = {CY\ and 
H G End{n). 

We will use the auxiliary space notation, where indices indicate on which sites (of the 
chain) the operators act non trivially. For instance, for an operator acting on two sites, 
O G End(C^ (X> C^), C34 will denote the operator O acting on sites 3 and 4 of the chain, i.e. 

O34 = I ® I ® C ® (I)®(^"^) G End(H) , (2.1) 

where I denotes the 2x2 identity matrix. 

Hence, in fll.ll) . B^ is the left boundary matrix B~^ acting on the first site 1: 

= (2.2) 

L-1 

hi^e-^-i are local Hamiltonians acting on sites {£, i + 1) 

hj^i+i = lO . . . . . . (8)1 , (2.3) 

and Bj^ is the right boundary matrix B~ acting on the last site L. 
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Pe/+i is the permutation operator between site i and site i + 1. As a matrix, the 
permutation P takes the form 



/ 1 








\ 








1 








1 








^0 








1 / 



(2.4) 



The Hamiltonian fll.ip describes the interaction of spins (up or down) among themselves, 
and with two boundaries describes by the matrices 5^. We are looking for its eigenstates 
and its spectrum. 

Contrary to the periodic case, the open XXX Hamiltonian does not possess an su{2) 
symmetry, so that the spin cannot be used as a quantum number. For non-diagonal boundary 
matrices, it does not commute with the u{l) generator, so that the pseudo-excitation number 
is not a good quantum number either. We will come back on this point in the next section. 



Reference state 

The state 

\t-.-t)e{cT' (2.5) 

is an eigenstate: 

i7|t...t) = (« + 7)lt...t). (2.6) 

This state will be chosen as a reference state (the so-called pseudo- vacuum) , and states with 
some spin down will be considered as pseudo-excitations above this reference state, see eq. 
(13. 2p . Let us stress that the reference state is not, in general, the ground state, and that 
spins down are not physical excitations. It is just a convenient way to parametrize all the 
states in "H. 

On contrary, due to the left boundary, the state | J, . . . |) is not an eigenstate. 



Hermiticity 

Note also that with this choice of boundaries, the Hamiltonian H is not Hermitian anymore 
for generic values of the boundary parameters. However, we will see that the energy is the 
same as the one computed for diagonal boundary matrices, so that the energy is real when 
the parameters a, /3, 7 and S are real. 

Remark that H can be pseudo-Hermitian: 

H^ = UHU^ with t/ = a^®...®a^ when a* = /3 ; (5* = 7 ; /i* = /i , (2.7) 

where * denotes complex conjugation. This case is still outside the range solved by usual 
coordinate Bethe ansatz (for 7^ 0). 
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Integrability 

The model we consider is integrable. It can be built from the following transfer matrix 



t u) 



- Tro(Ko^(n)To,<i...L>(n)i^o^(n)fo,<i...L>H), (2.8) 

To,<i...L>{u) = Rql{u) . . . Roi{u) with RQi{u) = Pqi + u , (2.9) 

To,<i...L>(n) = RM---Rol{u), (2.10) 

K-{u) = I + uB~ with B-=(^'~^ ^ ), (2.11) 



(5-7 

K+{u) = I + uB-^, with 5+=(^"~^ f3-a) ^^"^^^ 

where the K~{u) matrix (resp. K^{u) one) obey the reflection equation (resp. dual reflection 
one). Standard calculations |1] show that one recovers the Hamiltonian 



1 dt{u) 



i7 + (L-l-^±A±^)I. (2.13) 



u=0 



2 du 

3 Generalized Coordinate Bethe Ansatz 

The fact that the left boundary i?"*" is not diagonal anymore implies that the usual coordinate 
Bethe ansatz fails. Indeed, since B^ is triangular, it can flip a spin down to spin up, which 
is interpreted as the annihilation of a pseudo-excitation (or equivalently as its transmission 
outside the system). Hence, one cannot consider an eigenfunction with a given (fixed) number 
of pseudo-excitations. This is also another way to see that the spin is not a 'good' quantum 
number, since the Hamiltonian changes it. However, since there is only annihilation, one 
can consider eigenfunctions having a fixed maximum number of pseudo-excitations. It leads 
us to the following ansatz: 

n 

= E E E ^"'"^ 1^-+!' • • • '^") ' (3.1) 

■m=0 x',„+i< - <x„ geGm 

where Gm is a full set of representatives of the coset BCn/BCm {Gq = BGn, by convention) 
and BGm is the Bm Weyl group, generated by transpositions crj, j = 1, . . . , m — 1 that ex- 
change kj and /cj+i, and the reflection Ri exchanging ki and —ki. The vectors Ix^+i, . . . , x„) 
are given by 

|w,...,a;„) = |t...t i t-.-t i t t i t---t)G (C^)""^ (3.2) 

and we introduce the notation fc^™-* for the following truncated vector 

k^r^ = {kg{m+l),---,kg^ri)) ■ (3.3) 
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For this definition to be consistent, the coefficients Ag ' do not have to depend on the 
choice of the representative i.e. 

4r^ = 4"'"^ ioT any heBCm- (3.4) 

The coefficients Ag"'™'' are complex numbers to be determined such that $„ is an eigen- 
function of H i.e. such that the following equation holds 

H<^n = . (3.5) 

We project equation (13. 5p on the different independent vectors to get constraints on the 
coefficients aIP'"^\ 

Since if is a sum of operators acting on (at most) two neighboring sites only, one has to 
single out the cases where the x's obey the following constraints: 

• all the x/s are far away one from each other (1 + Xj < Xj+i, Vj) and are not on the 
boundary sites 1 and L. This case will be called genericH. 

• Xj + 1 = Xj+i for some j, 

• Xi = 1, or Xm = L. 

As the eigenvalue problem is a linear problem, it is enough to treat the cases where at most 
one of the particular cases appears: more complicated cases just appear as superposition of 
'simple' ones. 

Calculation of the energy: projection on \xi, . . . , x„) for (a;i, . . . , Xn) generic 

As in the usual coordinate Bethe ansatz [2], this projection provides the energy: 

E = a + + ^\{e'^^) where \{u) = n + -- 2= ^^~ ' . (3.6) 
j=i 

Let us remark that, up to the boundary terms a and 7, the energy takes the same form as 
in the periodic case. 

We want also to stress that the coefficient /i does not enter the energy, hence the spectrum 
for the model is the same as the spectrum of the model based on diagonal matricecl- 

^Here and below, unless explicitly specified, all the (sub)sets of XiS will be considered as generic. 
■^To be precise, one sfiould also check that the Bethe equations do not depend on yu cither: one will see 
below that it is indeed the case. 
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Scattering matrix: projection on \x\, . . . ,Xj, ^j+i = ^ -\- Xj, Xj^^i • • • ? Xn) 

As usual, this projection provides tlie scattering matrix between pseudo-excitations. It 
is given by a relation between Ag"'*^^ and Ag^f'' where cxj is the permutation of j and j + 1. 
Namely, we get 

A(^.o) = ^(e*^^'^ , e^*^«(^+i)) 4"'°) , (3.7) 

m X a(u,v) , , , .2v — uv — l ,^ . 

with b{u,v) = and a[u,v) = t . (3.8) 

a[v, u) uv — 1 

The normalization chosen for the function a(M, v) is for further simplifications. As expected, 
this relation is similar to the periodic case since the boundaries are not involved in this 
process. 

First relation: projection on |1, Xm+i • • • ? Xn), m > 1 
This relation is a new one. We get, for any g G G^, 

m 

^(n,m-i)^.fc^,,„ _ - 1 + e^^-«'^- - A(e*'=«0) = , (3.9) 

where Hm = BCm/ BCm~i- To obtain this relation, we have used the following property 
(valid for any function /): 

/(^)e^'=s(™)e*4'"'^''"' = Y e^'"^'"''^^'"' ^ f{gh)e'^^^<-^) . (3.10) 

Let us stress that equation (13. 9p does not depend on the choice of representative of Gm- 
Reflection coefficient for the left boundary 

The reflection coefficient is computed using relation (13. 9p at m = 1. In that case. Hi = 
BCi is constituted of 2 elements only: the identity id and the reflection Ri that changes ki 
into —ki. One gets 

45;'=fl(a"..)<'»' wth = -„^i_i±|Z^ = !im, (3,11) 

where r+(M) is given in ( 13.15^ . 

Remark that R{u) R{-) = 1: for the physical excitations, the left boundary is purely 
reflective (no loss of excitations). 
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Second relation: projection on |a3j„_|_i . . . , Xn), m > 1 



This projection provides a relation between the coefficients from the level m — 1 and m 
since we must take into account that the left boundary can destroy one pseudo-excitation 
present on the site 1. We obtain the following constraint, for any g G Gm, 

m 

4'^'""'^e'^«^(™) - J2 4"'"^ ^(e'^'O = . (3.12) 

h&Hm i=l 

Transmission coefficient for the left boundary 

From fl3.9p and (I3.12p . we may express all the Ap"'"*^ (m > 1) in terms of A^J^'^^ thanks 
to the recursive relation 

_ 2^(m)^gifc,i^ _ _ _ ^ ^ik,^^ j^{n,m-l) 1 < m < n , (3.13) 

with the following definitions: 



(m-l 
Y[a{um,Uj)a{uj,l/u^)\ , (3.14) 



l-u + (3-a (u- l){l-u + (3 - a) 

r+{u) = X{u) = . (3.15) 

1 — u[l + u) 

Relation f l3.13p can be interpreted as the transmission of one (among m) pseudo-excitation(s) 
through the left boundary: this pseudo-excitation has been destroyed (from the spin chain 
view point). 

Remark that T*^™) is proportional to /i, in accordance with the picture described at the 
beginning of the section: when /i = 0, the boundary becomes diagonal, the usual coordi- 
nate Bethe ansatz works, and there is no relation between A^P'"^^ and each level m 
providing an independent eigenfunction (with a fixed number of pseudo-excitations); when 
/i 7^ all the levels m are related, but we get an independent eigenfunction for each upper 
number n of pseudo-excitations. 

Iterative use of f l3.13p leads to 

j^in,m) ^ T(™)(e*'=«^...,e*'='"")4'''°^ (3.16) 

tM(«„...,«J = I n a(«.,«.)a(«„lM))(3.17) 



vi<i<^<m 



The proof that fl3.13H3.l4|) is a solution of both equations (13.90 and fl3.12p is postponed 
to section H] and relies on a residue computation. The integrability of the model plays a role 
at this place, since there are a priori too many constraints but not all of them are independent. 
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Note that relations (13 .7^ . f l3.1ip and ( I3.17P show that all the coefficients A^"'"^^ can be 
expressed in term of A^^'^^ . It is important to notice that this computation is consistent since 
the obtained expression does not depend on the way we write g in terms of the generators 
CTj and -Ri. 

Bethe equations: projection on \xi . . . , Xn-i, L) 

This last constraint consists in the quantization of the pseudo-excitations moments since the 
system is in a finite volume. In the context of the coordinate Bethe ansatz, this quantization 
leads to the so-called Bethe equations, explicitly given by 



e'"', e'"^) Sie-'^^^e"") = e''"-"^ 7_ ( ~; _./ for l<j<n (3.18) 



i=i 



rile \ r_(e 



{u) = - ^(l-n + (5-7). (3.19) 



u 



We remind that the scattering matrix S{u,v) is given in (13. 8p . while the r+{u) function is 
defined in (I3.15p . 

Let us remark that Bethe equations (13.181) do not depend on the parameter /i. Therefore, 
as previously mentioned, the spectrum (i.e. the energy (13. 6p ) is also independent of this 
parameter. Therefore, the spectrum is similar to the one with diagonal boundaries (yU = 0). 
In [25], similar results on the spectrum have been obtained via algebraic Bethe ansatz. In 
contrary, the eigenvectors depends on the parameter /i via relation (13.170 determining the 
coefficients entering in our ansatz. The eigenvectors computed thanks to the algebraic Bethe 
ansatz in [25] also depends on fi since the creation operators B{X), used to construct the 
ansatz, depends on /i. Unfortunately, a direct proof that the eigenvectors constructed via 
both methods are identical is a difficult task and beyond the scope of this paper. 



Completeness of the ansatz 

It is clear that two states $„i and are independent when ni 7^ Thus it remains to 
prove that the Bethe equations provide the right number of solutions, and that, at given n, 
these solutions are independent. 

Moreover, it is believed that the coordinate Bethe ansatz for open XXX chaiij^ is complete 
when fi = 0. Since the Bethe equations do not depend on fi, they provide the same number 
of solutions. Hence, from the conjecture at /i = 0, we deduce that when fi ^ 0, the set of 
solutions has the right dimension to get the complete set of eigenstates. 

"^To our knowledge it is only proven for closed spin chains, see e.g. |28 1 l29 j . 
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4 Proof of the transmission relation (13.131) 

In this section, we prove that ( 13.13^ imphes (13. 9 p and ( I3.12p . 

We start with f l3.12p . First, we remark that a consequence of f l3.13p is 



A(n,m) ^ A(n,m) I T^'") (e'''3i ,...,e''°g(™-l) ,e'''9("'+i) ) q/ 



1 < j < m — 1 , 

g«fcgm ^ g*'S5(m + l) ^ j = fjl ^ (4-1) 



9(i+i) 



j > m + 1 . 



Then, using again fl3.13p to express now in terms of A^J^'^' and using (14.1 p to express 



(n,m) 



^^^'^^ in terms of relation flXT^ becomes the functional relation 



I (n,m) 



m+1 



m+l 



m+1 



Ujr+{uj) J]^ a{uj,ue) a{ue, — ) + — r+(— ) J]^ a(— a{ue,Uj) - X{uj) 



£=1 



J J J i=\ J 



0,(4.2) 



where Uj stands for exp{ikgj) and the functions are defined in (13. 8p and (I3.15p . 
To prove this last relation (14.20 . let us introduce the following function 



u 



1 — u + (3 — a 
2(1 -m) 



m ^ 

|Ta(M,Mf) a{ui, -) . 



Looking at the poles of F^^\u)^ one can compute the residues of this function: 



m+1 



Res(F*-™'''(n)) = Uj r^{uj) TT a{uj, ui) a{ui, — ) 



=1 
m+1 



1 1 1 

Res(F(™)(n)) = - r+(-) TT a(n„ n,) a(n,, -) 



Res(F(™)(n)) 
Res(F(™)(^)) 



M=l 



--(/3-a) 



m+1 



^(/3-a)-^A(n,) 



(4.3) 

(4.4) 

(4.5) 

(4.6) 
(4.7) 



Then, (g^l) is equivalent to 'Eresidue F^^^H^) = which proves fl312D . 

From (13. 9p . we use the same procedure to obtain a new functional relation. After use of 
(14. 2p . this relation simplifies to 



m+1 



m+1 



m+1 



" 111 1 

J •'J.J. 7/. 7;^ 7/. J. J. 7/. 



^=1 



Uj U, Uj ^ ^ Uj 



rn rn 

+ (/^ - « - 1 - 5^ A(n,)) ( X{uk)) = . 

J=l k=l 



(4. 



Again, this relation is equivalent to a residue calculation, the function to consider being now 



5 Conclusion 

We have computed the spectrum and the eigenfunctions of the open XXX model with one 
triangular boundary matrix using a generalized coordinate Bethe ansatz. The next step 
would be to compute scalar products of the eigenstates to get informations on the correlation 
functions of this model. As far as scalar products are concerned, we think that an approach d 
la Gaudin [11] should work. A determinant formula, similar to the Slavnov determinant [30] 
will be needed to compute the correlation functions in a simple way. 

It remains to treat the case where the two boundary matrices are both triangular, and 
also the cases where one or two of the boundaries are general 2x2 matrices. Work is in 
progress on these two cases. We expect to find constrains between the boundary parameters 
in the first case. We remind that the method has been applied successfully in the case of 
open XXZ model [27]. The second case needs a further generalization of the coordinate 
Bethe ansatz. 

Finally, we believe also that the method presented here can be useful to find the spectrum 
for other integrable models with boundaries. 
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